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Abstract. 

The sunspot area fluctuations for the northern and the southern hemispheres of the Sun over 
the epoch of 12 cycles (12-23) are investigated. Because of the asymmetry of their probability 
distributions, the positive and the negative fluctuations are considered separately. The auto- 
correlation analysis of them shows three quasi-periodicities at 10, 17 and 23 solar rotations. The 
wavelets gives the 10-rotation quasi-periodicity. For the original and the negative fluctuations the 
correlation coefficient between the wavelet and the auto-correlation results is about 0.9 for 90% 
of the auto-correlation peaks. For the positive fluctuations it is also 0.9 for 70% of the peaks. For 
^^' 90% of cycles in both hemispheres the auto-correlation analysis of negative fluctuations shows 
that two longer periods can be represented as the multiple of the shortest period. For positive 
fluctuations such dependences are found for more than 50% of cases. 
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1. Introduction 



In the last decades, the intermediate quasi-periodicities of many solar activity tracers 
have been discussed. The about 12-rotation periodicity identified Krivova & Solanki 
(2002) for sunspot data during 1749-2001. It was prominent during times of stronger 
activity, whereas it diminished and sometimes faded into the background during weak 
cycles. Getko (2006) found it in both high and low activity periods for the monthly Wolf 
numbers during cycles 1-22 and for the group sunspot numbers during cycles 5-22. Two 
^y-N ] longer quasi-periodicities at 17 rotations and at 23 rotations were found in many solar 
^\ , activity parameters from the bottom of the convection zone to the atmosphere. More up- 
CO to-date review is by Obridko and Shelting (2007). Here I present results from a statistical 

l/~j study of these periodicities in the sunspot areas during cycles 12-23. It enables one to 

(^ ' deduce the mean length of the time period between strong fluctuations. 

o' 

'^ _ 2. A detailed analysis of quasi-periodicities 

I consider the daily sunspot areas for the northern hemisphere (Df), and the southern 
hemisphere (Df) for solar cycles 12-23 available at the National Geophysical Data Center 
C^ ■ (NGDC) (http://solarscience.msfc.nasa.gov/greenwch/). For the i-th Carrington rotation 
I evaluate the mean sunspot area for the northern hemisphere (S*"): 5" = j- X]!=i -^i"' 
where L is the number of days for the i-th rotation. I define the fiuctuation {F") of the 
mean sunspot area (S*") from the smoothed mean sunspot area: F" = S*" — S*," for i = 
l,...,N, where 5f = ^ Ej-tLe ^z'- Each of the time series {i^/^} and {Ff} contains 
N = 1706 elements. Both have almost the same probability distributions. Fig. [T^ shows 
the histogram of {F^} with a fitted Gaussian. The Kolmogorov-Lilliefors and Shapiro- 
Wilk tests reject the hypothesis of normality for them. Because each distribution has 
positive skew, the positive and the negative fiuctuations are considered separately. For 
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Figure 1. a Histogram of {F^} with a fitted Gaussian, b Top: Auto-correlation function (ct) 
of {-F"} for cycle 18. Middle: Same as for the upper curve, but for {-F"~}. Bottom: Same as 
for the upper curve, but for {-F""''}. The dotted lines represent two standard errors of each Ct 
function, c Dependence of T2 and rs on r for positive fluctuations. Lower solid curve represents 
the regression line for the points {t,T2) (dots). Upper solid curve represents the regression line 
for the points (tjTs) (triangles). Dashed lines represent the 95 per cent confidence interval for 
each regression line, d Same as for c, but for negative fluctuations. 



the northern hemisphere they can be defined as follows: 

where F^ ^ „„ _ J where 

FP where K" > ' ~ 1 FP where 
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It is known that for a time series which contains Gaussian white noise and a sinusoidal 
component a probability distribution is symmetric and the auto-correlation functions 
(ct) of that time series, of its positive fluctuations and of its negative fluctuations should 
be the same. The functions Cr of {F^}^ {F^^} f^^d {i^" } for one solar cycle are different 
(Fig. [TJj). The functions Cr of {F^} and {i^" } for cycle 18 have the significant global 
maxima at r = 11 rotations and smaller maxima at ri — k * t for fc = 2 and 3. In 54% 
of 24 cases (12 cycles in each hemisphere) the functions c^ of the original fluctuations 
have significant maxima for t G [7, 13]. In 30% cases the maxima for such r belong to the 
interval [Icr, 2cr]. For positive fluctuations this contribution is 50% and 46% respectively. 
For negative fluctuations in 92% cases the maxima at r e [7, 13] are significant. The 
mean value of all r G [7, 13] for which the maxima are significant is approximately 10 
rotations for all three fluctuation groups. I also consider the c,- maxima for r e [14, 19] 
and T e [20,27]. For more than 50% of cases the positive fluctuations create the auto- 
correlation peaks for which the periods ti « 17 and T2 k, 23 can be represented as 
Tk ^ k*T where r e [7, 13] and A: = 2 or 3. For each of k the points (r, t^), the regression 
line (solid) and the 95% confidence interval for each line (dotted) are shown in Fig. \Vp. 
For fc = 2 the correlation coefficient for 13 points is 0.91, for fc = 3 it is 0.86 for 8 points. 
For the negative fluctuations the strong dependence between the considered periodicities 
was found in ^ 90% for fc = 2 and in ^ 50% for fc = 3 (Fig.[TJi). For k — 2 the correlation 
is 0.95 for 23 points and for fc = 3 it is 0.91 for 11 points. It is important to add that the 
Cr values at r > 27 are not reliable because of the solar cycle length. 

I also applied the Morlet wavelet (Torrence & Compo 1998) to three fluctuation time 
series for each of 24 cases. Fig. [2] shows the normalized wavelet maps for cycle 18. Black 
contours denote the 95 per cent significance level for detected peaks. The wavelet map 
of {Fp} (top) shows two significant peaks (at r w 6 and 8). They are mainly created 
by three strong fluctuations at the begining and at the end of the high activity period. 
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Figure 2. a-c: Wavelet power spectra of a: {-F"}, b: {-F" } and c: {F"^} mapping a time-fre- 
quency evolution of about 10-rotation periodicity. Top values of wavelet power are denoted by 
gradual darkening. Black contours denote significance levels of 95 per cent for detected peaks. 
A cone of influence is marked by the dashed region, e-g: Corresponding global wavelet power 
spectra, d: Time series {-F"} for cycle 18. 

During the remaining part of this period a rised power is at r = 11 rotations. Moreover, 
the integrated spectrum (right) shows the maximum at r — 11 and confirms the auto- 
correlation results. For {F"^} (middle) the peak at r = 11 is well detected with 95 per 
cent level, extends in time during the high activity period and dominates the integrated 
spectrum. The map of j-F"^} (bottom) is similar to the map of {i^"}, but the peak at 
T = 11 is significant. The global spectrum shows two almost the same peaks at r = 11 
and 8. Such an analysis was done for all 24 cases. The auto-correlation and the wavelet 
results are similar for r G [7, 13] (the correlation between them is 0.9 for 87%, 92% and 
72% of the auto-correlation peaks of {F^}, {F"~} and {F^'^} respectively). 

These results could indicate that the 10-rotation quasi-period is dominant. Getko 
(2004) showed that large activity complexes were responsible for strong sunspot number 
fluctuations. Thus, the time between strong fluctuations of toroidal magnetic flux in the 
tachocline could be on the order of 7-13 rotations. It is also well known that two solar 
hemispheres show certain hemispheric asymmetries in their solar-cycle features. However, 
two-sample Kolmogorov-Smirnov test shows that the 10-rotation quasi-periods evaluated 
for each of 12 solar cycles in each hemispheres do not differ. Two longer quasi-periods 
at about 17 and 23 rotations could be treated as subharmonics of the 10-rotation quasi- 
period (Figs.[T]; and [Hi)- This facts could explain a wide range of periodicities in various 
solar indices at all levels from the tachocline to the Earth. 



3. Conclusions 

( a) For both hemispheres the probability distributions of fluctuations are similar and 
have an asymmetry which means that there are more negative than positive fluctuations. 

(&) The auto-correlation analysis of the original, the positive and the negative fluc- 
tuations prefers three quasi-periods: around 10, 17 and 23 rotations. The wavelet maps 
show one dominant quasi-period at about 10 rotations. 

(c) For 90% of solar cycles in both hemispheres the auto-correlation analysis of nega- 
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tive fluctuations gives peaks for which the period T2 ~ 17 rotations can be represented 
as Tfe « fc * T where r G [7, 13] and k = 2. For fc = 3 a such dependence are rehable in 
50% of considered cases. For positive fluctuations such dependences are found for more 
than 50% of solar cycles in each hemispheres. 
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